Charging induced asymmetry in molecular conductors 
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We investigate the origin of asymmetry in various measured current-voltage (I-V) characteristics 
of molecules with no inherent spatial asymmetry, with particular focus on a recent break junction 
measurement. We argue that such asymmetry arises due to unequal coupling with the contacts and 
a consequent difference in charging effects, which can only be captured in a self-consistent model 
for molecular conduction. The direction of the asymmetry depends on the sign of the majority 
carriers in the molecule. For conduction through highest occupied molecular orbitals (i.e. HOMO 
or p-type conduction), the current is smaller for positive voltage on the stronger contact, while 
for conduction through lowest unoccupied molecular orbitals (i.e. LUMO or n-type conduction), 
the sense of the asymmetry is reversed. Within an extended Hiickel description of the molecular 
chemistry and the contact microstructure (with two adjustable parameters, the position of the 
Fermi energy and the sulphur-gold bond length), an appropriate description of Poisson's equation, 
and a self-consistently coupled non-equilibrium Green's function (NEGF) description of transport, 
we achieve good agreement between theoretical and experimental FV characteristics, both in shape 
as well as overall magnitude. 



PACS numbers: 85.65.+h, 73.23.-b,31.15.Ar 

Future electronic devices are quite likely to incorpo- 
rate molecular components, motivated by their natural 
size, mechanical flexibility and chemical tunability. En- 
couraging progress in this direction has been achieved 
with the capacity to self-assemble, functionalize and re- 
producibly measure the current-voltage (FV) character- 
istics of small groups of molecules. Molecular FVs re- 
veal a wide range of conducting properties, from metallic 
conduction in carbon nanotubes [1] and quantum point 
contacts [2] , to semiconducting behavior in DNA [3] and 
conjugated aromatic thiols [4], and insulating behavior 
in alkylthiol chains [5] . Interesting device characteristics 
such as rectification [6] , switching [7] and negative differ- 
ential resistance on silicon substrate [8] have also been 
reported. In particular, molecular rectification continues 
to be a widely studied property, right from the earliest 
days of this field. 

The classic paradigm for asymmetry in molecular I- 
V measurements is the Aviram-Ratner diode, consist- 
ing of a semi-insulating molecular species bridging an 
electron donor-acceptor pair [9]. A positive bias on the 
contact at the donor end brings the energy levels on 
the donor and acceptor sites into resonance, while the 
opposite bias moves the system away from resonance, 
leading to a strongly asymmetric I-V characteristic [6]. 
Another example of asymmetric I-V is when the elec- 
trostatics of the system is dominated by one contact, 
e.g., a gated molecule where the gate is electrically con- 
nected to the source electrode [10,11]. This gives rise 
to a spatially asymmetric electrostatic potential over the 
molecule which in turn generates a strong asymmetry in 
the I-V. In these, as well as most commonly studied ex- 
amples of rectification, some kind of spatial asymmetry 



in the system seems essential, causing the energy levels, 
the electrostatic potential and the electron wave func- 
tions to be quite different for positive and negative volt- 
ages [12]. Typically, such asymmetry leads to peaks of 
similar heights in the conductance- voltage (G-V) charac- 
teristic occurring at different bias values. 

In this paper, we address an I-V asymmetry observed 
for a spatially symmetric molecule [13-15] that is quali- 
tatively different from and weaker than the rectifications 
in the above mentioned situations. The molecular I-V 
curves for these systems start off being symmetric, but 
pick up a weak, reversible asymmetry as the contacts 
are manipulated. In contrast to the two cases of recti- 
fication mentioned above, conduction in these molecules 
at opposite voltages occurs essentially through the same 
molecular levels with very similar wave functions. For 
resonant conduction in particular, this asymmetry shows 
up as conductance peaks of different heights occurring at 
symmetrically disposed voltage values. We show that the 
origin of the observed contact-induced asymmetry is non- 
trivial, and involves self-consistent shifts in the energy 
levels due to charging effects. Asymmetry in charging 
arises due to unequal coupling with the contacts, and 
seems to be present in conduction measurements per- 
formed with a break junction [13] or an STM tip [15]. 

Although this work addresses the origin of a weak rec- 
tification effect that may or may not be of practical sig- 
nificance from a device point of view, it serves three addi- 
tional purposes: (a) using a computationally inexpensive 
yet rigorous self-consistent transport model we show the 
significance of charging effects in explaining experimen- 
tally observed I-V characteristics; (b) our calculated I-V 
(Fig. 1) is not just qualitatively, but also quantitatively 
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in agreement with break-junction measurements on very 
small groups of molecules [13]. This is a significant devi- 
ation from the typical situations where theoretical esti- 
mates of current values differ from experimental currents 
by orders of magnitude [16-19]; (c) the physics also sheds 
light on the nature of the charge carriers, i.e. whether 
conduction is n-type or p-type. The identification of the 
polarity of charge carriers is of obvious importance in 
semiconductor devices, since electrons and holes have dif- 
ferent effective masses, leading to different mobilities. By 
analogy, molecular LUMO and HOMO levels have quite 
different wavefunctions [20], leading to different trans- 
missions and current conducting properties. Much un- 
certainty exists about the nature of the conduction or- 
bitals, or equivalently, the position of the Fermi energy 
relative to the molecular energy levels [17,16,19,21,22]. 
We argue that for the kind of contact-induced asymme- 
try discussed here, current is lower for positive bias on 
the strongly coupled contact for HOMO conduction, and 
higher if conduction is through a LUMO level. 

A description of molecular conduction can be put on 
a quantitative footing by using an appropriate molecular 
Fock matrix F (ab-initio or semi-empirical), coupled with 
a non-equilibrium Green's function (NEGF) formulation 
of transport [17,23]. For a given Fock matrix F, overlap 
matrix S, and contact self-energies Si. 2 with correspond- 
ing broadenings Ti t 2 = i (^1,2 — ^1.2) > tne energy lev- 
els are given by the poles of the noncquilibrium Green's 
function G, while their occupancies are obtained from 
the corresponding density matrix p, the contact electro- 
chemical potentials p\ 2 and the Fermi functions f\$: 

G(E) = (ES - F -£ a -£ 2 ) _1 

/OO 
dE (/iGr 1 G t + / 2 GT 2 G t ) 
-co 

fiAE) = [1 + exp {{E - Mi^AsT)]- 1 (1) 

The number of electrons N and the steady-state current 
I are then given by: 

N = 2 (for spin) x trace(p S) 
2e r°° 

I=- dET{E,V)[h (E) - f 2 (E)] (2) 
h J-00 

where the transmission T is given by: 

T(E, V) = trace(riGr 2 G t ) (3) 

It has been pointed out that the nature of the molecu- 
lar I-V could sensitively depend on the self-consistent po- 
tential profile [24]. In our calculations, the self-consistent 
charging effect is included within an extended Hiickel 
(EHT) [24,20] description of the Fock matrix F, coupled 
with the NEGF equations described above. The self en- 
ergies are calculated for Au (111) contacts in EHT using 
a recursive technique [20]. Both charging and screening 



effects are incorporated in F through a self-consistent po- 
tential Use = y{f>) which describes the potential profile 
of the molecule under applied bias: 

V(A P ) = V Lap i ace + AV P 

oisson 

(Ap) + AV mage (Ap) (4) 

where Ap represents the change in density matrix under 
bias, p — peg- Equations (1) and (4) are then solved self- 
consistently to obtain the converged Fock matrix F. Any 
asymmetry in the potential profile is included explicitly 
in the self-consistent potential Use- 
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FIG. 1. (Color online) I-V characteristics of the 
gold-molecule-gold system shown above. Solid line: theo- 
retical calculations; dots: experimentally obtained data in a 
break junction configuration [13]. The molecular energy levels 
are raised by a constant potential V c (treated as a fitting pa- 
rameter) relative to the gold Fermi energy (-9.5 eV) in order 
to simulate the effect of charge transfer and band line-up at 
equilibrium. The same value of V c (1.5 eV) is used to gener- 
ate all three curves, although allowing slight variation gives a 
better fit. The upper curve (a) is obtained on decreasing the 
left electrode coupling by stretching the sulphur-gold bond 
length from its equilibrium value of 2.53 A to 3.18 A and the 
lower curve (c) corresponds to the reverse, whereas curve (b) 
represents the symmetric coupling situation. The bias polar- 
ity is defined as positive when the applied voltage on the left 
contact is positive. Curves (a) and (c) are offset by +I/1A 
and -lpA respectively for better visibility. 

The Poisson part in equation (4) is approximated us- 
ing the CNDO (complete neglect of differential overlap) 
method [25,26], of which only the Hartree potential is be- 
ing utilized in our treatment. Both charging and screen- 
ing effects are incorporated into this term. In CNDO 
approximation, this Poisson term becomes: 

AlC isson (Ap) = (A p )m7m + ( A p)a7ma 

A 
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where (Ap) 



M 



M 



A I 



7m = (mmImm) 

1m a = (mmIMa) 
M, A = Atomic sites 

(i = Slater type atomic orbital 



(5) 



The notations used in Equation (5) are consistent with 
[26]. The two electron integrals 7's, are the CNDO pa- 
rameters which are obtained from the experimental data 
and empirical fitting. These CNDO parameters allow us 
to capture the main physical characteristics of the e-e 
interaction without evaluating the computationally ex- 
pensive two electron integrals directly. 
,0 
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FIG. 2. (Color online) Equilibrium transmission coefficient 
as a function of energy at equilibrium. Solid line: symmet- 
ric coupling; dotted line: asymmetric coupling with the left 
contact weakly coupled. The transmission coefficients do not 
show any significant difference in the energy range of interest 
from -10 eV to -9 eV. As current is proportional to transmis- 
sion, this explains why current values do not differ appreciably 
although the S-Au bond length has been stretched by 0.65 A 
in the asymmetric coupling situation. 

The Laplace and image potentials are calculated by us- 
ing a finite element method [27] treating the atomic sites 
as points in free space. The two electrodes are treated as 
metallic plates [100 A x 100 A] separated by the molec- 
ular length. The molecule is placed in between the two 
plates. The Laplace part is then obtained by solving 
Laplace's equation in 3D with the boundary conditions 
set by jjL\ t 2 — Ef=f0.5V ap pi (V app i: applied bias). The im- 
age part is calculated similarly but with different bound- 
ary conditions set by the potentials on the metallic plates 
due to point charges on the atomic sites. The inclusion of 
image potential in our model does not make any signifi- 
cant differences other than lowering the charging energy 
of the system by around 0.1 eV. 

To illustrate charging induced I-V asymmetry we ap- 
ply our self-consistent transport model to the sym- 



metric molecule, [9,10-Bis((2 -para-mercaptophenyl)- 
ethinyl)-anthracene] (see top of Fig. 1) studied in a re- 
cent break junction measurement [13]. The curve (b) in 
Fig. 1 is obtained assuming ideal sulphur-gold bonding 
(2.53 A) on both sides of the molecule (symmetric cou- 
pling situation) on a gold-molecule-gold system. On the 
other hand, the curves (a) and (c) represent an asymmet- 
ric coupling situation where the asymmetry is introduced 
in the system by stretching the bond between sulphur and 
gold on the left side and the right side respectively by 0.65 
A. The sulphur-gold bond length has been treated as a 
fitting parameter which is justifiable as the exact bond 
length is not known from the experiment. 




(a)V>0 



(b)V<0 



FIG. 3. (Color online) Origin of asymmetry due to charg- 
ing shown for an one level model: one side (say the right sub- 
strate) is strongly coupled. Although the same level is crossed 
by the contact electrochemical potentials Qtti and ^2) in both 
bias directions, for negative applied bias (b) the HOMO level 
is emptied out by the strongly coupled contact, which posi- 
tively charges the molecule and shifts the energy levels down. 
Such a shift, not present for positive applied bias (a), post- 
pones the onset of conduction and effectively stretches out the 
voltage axis in the I-V along the direction of negative applied 
bias. 

Another fitting parameter used in our calculation is a 
constant potential V c which is applied on the molecule 
(only on the whole molecule and not on any part of the 
gold contacts) in order to move the molecular energy 
levels up or down rigidly relative to the contact Fermi 
energy. The precise location of equilibrium Fermi en- 
ergy Ef can depend sensitively on many factors such as 
surface conditions, environmental situations, geometrical 
fluctuations etc. In the absence of detailed information 
of all these factors in a particular experimental set-up, it 
is necessary to treat the position of Ef as an adjustable 
parameter. Thus we believe it is justifiable to include this 
adjustable parameter V c in our calculations to take into 
account of the effect of chrage transfer at equilibriumin 
and adjust the position of Ef relative to the molecular 
energy levels. Our fitting requires that Ef be closer to 
HOMO rather than LUMO i.e. the conduction is p-type. 
The Fermi energy of gold (within the EHT description) 
is kept fixed at -9.5 eV [20] and relative to this value the 
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molecular energy levels are moved up by V c — 1.5eV so 
that Ej — Ehomo is se t to be 0.33 eV. 
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FIG. 4. (Color online) (a) Change in the number of elec- 
trons and (b) the molecular energy levels as a function of ap- 
plied bias. Solid line: symmetric coupling; dotted line: asym- 
metric coupling with the left contact weakly coupled. Both 
the number of electrons and the energy levels change symmet- 
rically with applied bias in the symmetric coupling situation, 
as expected. On the other hand, for the asymmetric coupling 
case more electrons empty out in the negative bias direction 
(a) making the charging asymmetric and moving down the 
energy levels more in the corresponding bias direction (b). 

The calculated trend in the I-V characteristics agrees 
both in shape and current magnitude with the experi- 
mental result, as evident from Fig. 1. The agreement 
of our calculated current value with experiment suggests 
that the observed FV is obtained for a near ideal gold- 
molecule-gold contact and the current is going through a 
single or at most a few number of molecules. In spite of 
stretching the S-Au bond by 0.65 A the amplitude of cur- 
rent does not change appreciably. This can be explained 
by noting that the transmission coefficients do not dif- 
fer significantly in the energy range of interest from -10 
eV to -9 eV (see Fig. 2). The little discrepancy in the 
shape of the I-V can be attributed to the variation in 
the experimental results. It should also be noted that 
the same value for the constant potential V c (i.e. same 
Ef — Ehomo) is used to generate all the three curves in 



Fig. 1, though the position of the Fermi energy can con- 
ceivably be different for curves (a) and (c) as the experi- 
mental conditions are different. Indeed, by using slightly 
different values for V c within a justifiable range ( V c = 
1.55 eV and V c = 1.6 eV for curves (a) and (c) respec- 
tively) we do obtain a better fit of the I-V (not shown) . 

In our calculations we assumed that the system is in 
the self-consistent field regime. This assumption should 
be valid as long as the energy level broadening due to 
contact coupling is comparable to (or greater than) the 
single-electron charging energy. We estimate these fac- 
tors to be of the order of 0.2 eV and 1.5 eV, respectively. 
It should be further noted that the experimental current 
value is higher (of the order of /xA) than other typical ex- 
perimental values, suggesting strong chemisorbed bond- 
ing between sulphur and gold at both ends. Thus we 
believe our assumption of self-consistent field approach 
is a correct one. 



a) Symmetric b) Asymmetric 

coupling coupling 
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FIG. 5. (Color online) A color plot of the transmission as 
a function of energy and applied bias. Fig. 5a corresponds 
to symmetric coupling situation (Fig. lb) whereas Fig. 5b 
corresponds to the asymmetric coupling situation where the 
left side is weakly coupled (Fig. la). The white dotted lines 
pinpoint the onset of level crossing in both bias directions 
(enlarged portions shown below). The Fermi energy of the 
device (-9.5 eV) is assumed to be closer to the HOMO level 
and only this level is conducting in the applied bias range. 

In order to explain the origin of this asymmetry in 
the I-V we need to understand the effect of asymmet- 
ric charging on molecular conduction. This is explained 
schematically in Fig. 3 for a simple single level model. 
Let us suppose that the contacts are strongly asymmetric 
Ti). For positive bias on the weak contact (Fig. 
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3a), the latter is trying to empty the nearest (HOMO) 
level, while the strong contact is trying to fill it, with 
the net result that the HOMO level stays filled. Cur- 
rent onset occurs at the voltage where fi\ first crosses 
the (neutral) molecular level. For opposite bias direc- 
tion (Fig. 3b), however, the HOMO level is emptied 
out, which charges up the molecule positively. This adds 
a self-consistent charging energy that lowers the energy 
level, postponing thereby the point where the HOMO is 
crossed by /i2- In effect, this stretches out the voltage 
axis in the FV, leading thus to a smaller conductance for 
negative bias on the weak contact. For LUMO-based con- 
duction the argument is reversed, since filling the LUMO 
level charges up the molecule negatively. 

It can be seen from Fig. 4 that the number of elec- 
trons and the energy levels of the molecule change sym- 
metrically with the applied bias for the symmetric cou- 
pling situation. This is expected for a symmetric sys- 
tem. However, in the asymmetric coupling situation the 
weakly coupled left contact can not fill up the HOMO 
level quickly enough and so this level remains more empty 
in the negative bias direction compared to that in the 
positive bias (Fig. 4a). This in turn makes the charging 
asymmetric and forces the energy levels to move down 
more in the negative bias direction (Fig. 4b). 

Fig. 5 shows a color plot of the transmission T(E,V), 
for symmetric coupling (a) and asymmetric coupling (b), 
corresponding respectively to curves (b) and (a) in Fig. 
1. In both cases, the same HOMO level is crossed by 
the contact electrochemical potentials fi\ t 2 (white solid 
lines) leading to current conduction. For symmetrically 
coupled contacts, current starts flowing around ±0.8 V 
in both bias directions, leading to the symmetric I-V in 
Fig. lb. For the asymmetrically coupled case, however, 
differential charging leads to a postponed onset of con- 
duction at around -0.9 V for negative bias on the weak 
(left) contact, and an earlier onset at around +0.7 V for 
opposite bias polarity. This makes the current lower and 
stretched out in the negative bias direction which leads 
to the asymmetric I-V in Fig. la. The delay in the on- 
set is seen clearer in the enlarged portions shown at the 
bottom of Fig. 5. 

The charging induced asymmetry in the I-V also sheds 
light on whether the current conduction is p-type or n- 
type. If we would assume LUMO conduction in our cal- 
culations, the direction of asymmetry in the I-V would 
get reversed. Thus, from the knowledge of the cou- 
pling asymmetry and the I-V asymmetry we can predict 
whether it is HOMO or LUMO conduction. Although 
we do not know which contact is weaker for the break 
junction measurements, for similar measurements with 
STM [15] where the weaker contact is easily identified, 
this physics helps infer that conduction is p-type. Most 
often it is very difficult to pinpoint the position of Fermi 
energy relative to the molecular energy levels and so this 
ability of predicting conduction type is very desirable and 



useful. It should be noted that the nature of conduction 
could possibly be identified by other techniques, e.g. by 
using thermoelectric [28] or gating effects [29]. 

It is worthwhile pointing out that the I-V asymmetry 
discussed in this paper can not be explained solely in 
terms of asymmetry in the Laplace potential, neglecting 
charging effects altogether. Asymmetric charging gives 
conductance peaks of different heights at same voltage 
values (see Fig. 10 in [15]), in contrast to Laplace poten- 
tial asymmetries which give conductance peaks at differ- 
ent voltages but of the same height. 

In summary, we have achieved a quantitative agree- 
ment with experimentally observed asymmetric I-Vs and 
established that even for spatially symmetric molecules, 
contact asymmetry can induce an asymmetric I-V 
through asymmetric charging. The sense of the asymme- 
try depends on whether conduction is through a HOMO 
or a LUMO level. 
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